The  Laplace  transform  is  frequently  encountered  in  mathematics,  physics,  engineering  and 
other  areas.  However,  the  spectral  properties  of  the  Laplace  transform  tend  to  complicate 
its  numerical  treatment;  therefore,  the  closely  related  “truncated”  Laplace  transforms  are 
often  used  in  applications.  In  this  dissertation,  we  construct  efficient  algorithms  for  the  eval¬ 
uation  of  the  singular  value  decomposition  (SVD)  of  such  operators.  The  approach  of  this 
dissertation  is  somewhat  similar  to  that  introduced  by  Slepian  et  al.  for  the  construction 
of  prolate  spheroidal  wavefunctions  in  their  classical  study  of  the  truncated  Fourier  trans¬ 
form.  The  resulting  algorithms  are  applicable  to  all  environments  likely  to  be  encountered 
in  applications,  including  the  evaluation  of  singular  functions  corresponding  to  extremely 
small  singular  values  (e.g.  lO-1000). 
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Chapter  1 


Introduction 


The  Laplace  transform  £  is  a  linear  mapping  L2( 0,  oo)  —>  L2{ 0,  oo);  for  a  function  /  G  L2( 0,  oo), 
it  is  defined  by  the  formula: 

(A/))  (w)  =  J  e~tu]f(t)dt.  (1.1) 

As  is  well-known,  C  has  a  continuous  spectrum,  and  £_1  is  not  continuous  (see,  for  example, 
[lj).  These  and  related  properties  tend  to  complicate  the  numerical  treatment  of  C. 

In  addressing  these  problems,  we  find  it  useful  to  draw  an  analogy  between  the  numerical 
treatment  of  the  Laplace  transform,  and  the  numerical  treatment  of  the  Fourier  transform 
for  a  function  /  G  L1(M),  the  later  is  defined  by  the  formula: 

(A/))  M  =  (1-2) 

where  w£l. 

In  various  applications  in  mathematics  and  engineering,  it  is  useful  to  define  the  “truncated” 
Fourier  transform  Tc  :  L2(—  1, 1)  — >  L2(— 1, 1);  for  a  given  c  >  0,  Tc  of  a  function  /  G  L2(— 1, 1) 
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is  defined  by  the  formula: 


faf))(u)  =  j\~ict“f(t)dt.  (1.3) 

The  operator  Tc  has  been  analyzed  extensively;  one  of  the  most  notable  discoveries,  made  by 
Slepian  et  al.  in  1960,  was  that  the  integral  operator  Tc  commutes  with  a  second  order  differ¬ 
ential  operator  (see  El)  This  property  of  J-c  was  used  in  analytical  and  numerical  investigation 
of  the  eigendecomposition  of  this  operator,  for  example  in  [3]  and 

For  0  <  a  <  b  <  oo,  the  linear  mapping  Ca ^  :  L2(a,  b )  — >  L2{ 0,  oo),  defined  by  the  formula 

(£a,b(f))  (w)  =  [  e~tuJ  f(t)dt,  (1.4) 

J  a 

will  be  referred  to  as  the  truncated  Laplace  transform  of  /;  obviously,  Ca,b  is  a  bounded  compact 
operator  (see,  for  example,  CD- 

Bertero  and  Griinbaum  discovered  that  each  of  the  symmetric  operators  (£a,b)*  o  Lah  and 
£a,b  o  (£a,b)*  commutes  with  a  differential  operator  (see  [5]).  These  properties  were  used  in  the 
analysis  of  the  truncated  Laplace  transform  (see  0,  0)- 

Despite  the  result  in  [5] ,  more  is  known  about  the  numerical  and  analytical  properties  of 
Tc  than  about  the  properties  of  Ca,b- 

In  this  dissertation,  we  introduce  an  algorithm  for  the  efficient  evaluation  of  the  singular 
value  decomposition  (SVD)  of  £a,b,  and  analyze  some  of  its  properties.  A  more  detailed  analysis 
of  the  asymptotic  properties  of  Ca,b  will  be  presented  in  a  separate  paper. 

The  dissertation  is  organized  as  follows.  Chapter  [2]  summarizes  various  standard  mathe¬ 
matical  facts  and  certain  simple  derivations  that  are  used  later  in  this  dissertation.  Chapter 
[2]  also  contains  a  definition  of  the  SVD  of  the  truncated  Laplace  transform  and  a  summary  of 
some  known  properties  of  the  truncated  Laplace  transform.  Chapter  [3]  contains  the  derivation 
of  various  properties  of  the  truncated  Laplace  transform,  which  are  used  in  the  algorithms. 
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Chapter  [4]  describes  the  algorithms  for  the  evaluation  of  the  singular  functions,  singular  values 
and  associated  eigenvalues.  Chapter  [5] contains  numerical  results  obtained  using  the  algorithms. 
Chapter  [6]  contains  generalizations  and  conclusions. 


Remark  1.1.  Some  authors  define  the  truncated  Laplace  transform  as 
a  =  0,  or  define  the  operator  as  a  linear  mapping  L2(a,  b)  ->  L2(a,  b ).  See, 


in  (1.4),  but  allow 


for  example,  [7]. 
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Chapter  2 


Mathematical  preliminaries 


In  this  chapter  we  introduce  notation  and  summarize  standard  mathematical  facts  which  we 
use  in  this  dissertation.  In  addition,  we  present  a  brief  derivation  of  some  useful  facts  which 
we  have  failed  to  find  in  the  literature. 


2.1  Legendre  Polynomials 

Definition  2.1.  The  Legendre  polynomial  P of  degree  k  >  0,  is  defined  by  the  formula 


Pk(x) 


1  /  2  _  id 

2 kk\  dxk  X 


(2.1) 


As  is  well-known,  the  Legendre  Polynomials  of  degrees  k  =  0, 1....  form  an  orthogonal  basis 
in  L2(  —  1, 1).  The  following  well-known  properties  of  the  Legendre  polynomials  can  be  found 
inter  alia  in  0,  0: 


SjPk{x))2Ax  =  WTi 


(2.2) 
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(k  +  l)Pk+\(x)  =  (2k  +  1  )xPk(x)  -  kPk-i(x) 


(2.3) 


(1  -  x2)-^-Pk(x)  =  -kxPk(x)  +  kPk-i(x)  (2.4) 

dx 

^(1  -  x2)-^Pk(x)\  =  ~k(l  +  k)Pk(x)  (2.5) 

(2k  +  1  )Pk(x)  =  (Pk+i(x)  -  Pk-i(x))  (2.6) 

Po(x)  =  1  (2-7) 

P\(x)  =  x  (2-8) 

For  all  k>\1 


(k  +  l)Pfc+i(x)  =  (2k  +  1  )xPk(x)  -  kPk_i(x)  (2.9) 

In  this  dissertation  we  will  analyze  functions  in  L2(0,1);  it  is  therefore  convenient  to  use 
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the  shifted  Legendre  polynomials,  which  are  defined  on  the  interval  (0, 1). 


Definition  2.2.  The  shifted  Legendre  polynomial  of  degree  k  >  0,  which  we  will  be  denoting 
by  Pf,  is  defined  via  the  Legendre  polynomial  Pk  by  the  formula 


P*k(x)  =  Pk(  2x-l). 


(2.10) 


Clearly,  the  polynomials  Pf  form  an  orthogonal  basis  in  L2(0, 1).  The  following  properties 
of  the  shifted  Legendre  polynomials  are  easily  derived  from  the  properties  of  the  Legendre 


polynomials  by  substituting  (2.10)  into  (2.2  2.7). 


f„  (p‘(l))2dl  =  2Fn 


(2.11) 


xPk(x) 


1  (kPU(x)  i  d* /  \  i  (l  +  k)P*k+1(x)\ 

2fTT2 V+P*ix)  + - 1+2  k  ) 


(2.12) 


x(l 


k{  1  +  k) 
2(1  +  2  k) 


{P*k- rW  -  Pk+ 1(*)) 


(2.13) 


_d_ 

dx 


x) 


-k(l  +  k)P£(x) 


(2.14) 


Po(x)  =  1 


(2.15) 
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As  is  evident  from  (2.2)  and  (2.11),  neither  the  Legendre  polynomials  nor  the  shifted  Leg¬ 


endre  polynomials  are  normalized.  In  the  discussion  of  the  space  of  functions  L2( 0, 1),  we  will 
find  it  convenient  to  use  the  orthonormal  basis  of  the  functions  Pf(x). 

Definition  2.3.  We  define  P£ (x)  by  the  formula: 

P*(x)  =  P%{x)V2k  +  l,  (2.16) 


where  k  =  0,1, . 

Clearly,  the  polynomials  Pf  are  an  orthonormal  basis  in  L2(0, 1). 

Observation  2.4.  Pg  is  a  constant 

P*(x)  =  1  (2.17) 


Observation  2.5.  The  derivative  of  Pf  is  a  linear  combination  of  Pf ,  where  l  <  k.  The 


following  expressions  for  the  derivative  are  easily  verified  using  (2.6),  (2.16)  and  (2.10)  : 


d 

dx 


_ 3- 1  _ 

if  (x)  =  2v'2(2j)  +  l£  VWl  +  1)  +  1 


«=o 


(x) 


(2.18) 


i  J  x 

~^P2j+i{x)  =  2\/2(2j  +  1)  +  1  \/2(2l)  +  1  P2i(x ) 


(2.19) 


2.2  Legendre  Functions  of  the  second  kind 

Definition  2.6.  The  Legendre  function  of  the  second  kind  Qk{z)  is  defined  by  the  formula 

1  71 

Qk(z)  =  -  J  {z  -  t)_1Pfc(i)dt,  (2.20) 
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where  Pk(t)  is  defined  in  (2.1). 


The  following  identities  can  be  found,  for  example,  in  0,  0: 


Qk(z)  =  (-1  )k+1Qk(-z), 


(2.21) 


fc 

Qk{z)  =  / 
Jo 


d  4> 


\  fc+i ' 


(^z  +  \J  z1  —  1  cosh  (0)^ 


(2.22) 


Having  defined  the  shifted  Legendre  polynomials,  we  find  it  convenient  to  also  define  a 
shifted  version  of  the  Legendre  function  of  the  second  kind. 


Definition  2.7.  We  define  the  shifted  Legendre  function  of  the  second  kind  of  degree  k,  which 
we  will  be  denoting  by  Q*k,  by  the  formula 


Q*k(z)  =  Qk(2z  -  l) 


By  (2.16),  (2.20),  (2.21)  and  (2.23), 


(x  +  y)~1P*(x)dx  =  2(-l)kQ*k(y  +  l)V2k  +  l  y>  0 


(2.23) 


(2.24) 


and 


Q\(l  +  5/2)  =  Q\{l  +  5) 


f 


d 


\  fc+i 


^(1  +  5)  +  a/(1  +  S)2  -  1  cosh(0)^ 


(2.25) 


For  a  given  x  >  1,  Q%(x)  decays  rapidly  as  k  grows.  The  following  lemma  gives  an  upper 
bound  for  \Q*k(z)\,  where  z  >  x,  as  k  grows. 


Lemma  2.8.  Let  6  >  0.  We  introduce  the  notation  6  =  ^/(l  +  <5)2  —  1.  Then,  for  all  y  >  0 


\Q%0  +  +  y)\  < 


1 


(2.26) 


where  Q*k  is  defined  in  (2.23). 


Proof.  By  (2.25), 


\Q*  kO-  +  +  y)  |  — |Qfc(l  +  <5  +  2y)|  — 

ro 

Jo 


d  <)> 


2  —  1  cosh 


fc+i ' 


(2.27) 


Since  (1  +  h  +  y)  >  (1  +  5), 


I Q*kO  +  +  y)l  — IQfc(i  +  <5  +  2y)|  < 


< 


f 


d  (j) 


(2.28) 


2  —  1  cosh 


fc+1  ' 


Clearly,  5  =  y/(l  +  h)2  —  1  >  0,  and  by  (2.22) 


\Q*kO  +  <^/2  +  y)\  <  f 

Jo 


df 


0  ( 1  +  <5cosh(</>)V 


(2.29) 


We  define 


v  =  log 


(2.30) 
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and  break  the  integral  in  (2.29)  into  integrals  on  the  two  intervals  [0,  v)  and  [u,  oo): 


\Q*kO-  +  +  y)\  < 


dc/) 


0  (l  +  <5  cosh 


fc+i 


+ 


r  oo 

Ju  ( 1 


d0 


+  5  cosh 


fc+i ' 


(2.31) 


Clearly, 


1  +  6  cosh 


\  fc+i  — 


(*)) 


l  +  <5 


fe+i ' 


(2.32) 


and 


1 


< 


1 


\ k+l  —  /  \  ^ I ^ ^ 

1  +  <5  cosh(^)J  ( <5exp((/>)/2j 


(2.33) 


so  that, 


IQ*fc(l  +  <V2  +  y)|  < 


(2.34) 


Substituting  (2.30)  into  the  last  inequality,  we  obtain 


IQfc(l  +  <^/2  +  y)\  < 


1 


+ 


and  from  it,  we  obtain  (2.26). 


(2.35) 


a 
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2.3  Laguerre  functions 


Definition  2.9.  The  generalized  Laguerre  polynomial  (x)  of  order  a  >  —  1  and  degree 
k  >  0,  is  defined  by  the  formula 


(2.36) 


Definition  2.10.  The  Laguerre  polynomial  Lk(x)  is  the  generalized  Laguerre  polynomial  of 
order  0: 


Lk(x)  =  Lf\x) 


(2.37) 


As  is  well-known,  the  Laguerre  polynomials  are  an  orthonormal  basis  in  the  Hilbert  space 
induced  by  the  inner  product 


poo 

( f,g)=  /  e~xf(x)g(x)da 
■Jo 


(2.38) 


The  following  well-known  properties  of  the  generalized  Laguerre  polynomials  can  be  found, 
inter  alia,  in  [8j: 


Lk  1(x)  =  Lk(x)-Lk- 1(*) 


(2.39) 


d_ 

dx 


Lk(x)  =  -42  r 


(2.40) 
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xLk(x )  =  —(k  +  l)Lk+i(x)  +  (2k  +  1  )Lk(x)  -  kLk_i(x) 


(2.41) 


e  xtLk(x)dx 


(t  -  1  ft^-1 


(2.42) 


Lk( 0)  =  1 


(2.43) 


L0(x)  =  1 


(2.44) 


L\(x)  =  1  —  x  (2-45) 

For  all  k  >  1, 

(fc  +  l)Lk+i(x)  =  (2k  +  1  -  x)Lk(x)  -  kLk_i(x )  (2-46) 

It  is  convenient  to  use  functions  which  are  orthonormal  in  the  standard  L2(0,oo)  sense. 
Therefore,  we  will  use  the  Laguerre  functions,  as  defined  below,  rather  than  the  Laguerre 
polynomials. 

Definition  2.11.  We  define  the  Laguerre  function,  which  we  will  be  denoting  by  &k,  via  the 
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formula 


$k(x)  =  e  x/2Lk(x). 


(2.47) 


Clearly,  the  Laguerre  functions  &k(x)  are  an  orthonormal  basis  in  the  standard  L2(0,  oo) 


sense. 


Observation  2.12.  The  derivative  of  a  Laguerre  function  of  degree  A;  is  a  linear  combination 
of  Laguerre  functions  of  degree  k  and  lower.  The  following  expression  is  easy  to  verify  using 


(2.40)  and  (2.47): 


x  i= o 


(2.48) 


2.4  The  complete  elliptic  integral 

Several  slightly  different  definitions  of  the  complete  elliptic  integral  of  the  first  kind  can  be 
found  in  the  literature.  In  this  dissertation,  we  will  use  the  following  definition. 

Definition  2.13.  The  complete  elliptic  integral  of  the  first  kind,  K{m )  is  defined  by  the  formula 


K(m) 


—  msin2(0))  1//2  d 9. 


(2.49) 


2.5  Singular  value  decomposition  (SVD)  of  integral  operators 

The  SVD  of  integral  operators  and  its  key  properties  are  summarized  in  the  following  theorem, 
which  can  be  found  in  m- 

Theorem  2.14.  Suppose  that  the  function  K  :  (c,d)  x  (a,  b)  — >  M  is  square  integrable,  and  let 
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T  :  L2{a,b)  ->  L2(c,d)  be 


(T(f))(x)=  [  K(x,t)f(t)dt.  (2.50) 

J  a 

Then,  there  exist  two  orthonormal  sequences  of  functions  un  :  (a,  b)  — >  M  and  vn  :  (c,  d)  — >  M 
and  a  sequence  sn  E  M,  /or  n  =  0,  ...oo,  such  that 

OO 

K(x,t)  =  y ^Vn(x)snun(t)  (2-51) 

n=0 

and  that  so  >  si  >  ...  >  0.  The  sequence  sn  is  uniquely  determined  by  K .  Furthermore,  the 
functions  un  are  eigenfunctions  of  the  operator  T*  o  T  and  the  values  sn  are  the  square  roots 
of  the  eigenvalues  of  T*  o  T. 

Observation  2.15.  The  function  K  can  be  approximated  by  discarding  of  small  singular  values 
(see  PH): 

p 

K(x,t)  ~  Y.Vn^SnUnit)  (2.52) 

n=0 

2.6  Tridiagonal  and  five-diagonal  matrices 

In  this  section,  we  briefly  describe  a  standard  method  for  calculating  eigenvectors  and  eigen¬ 
values  of  symmetric  tridiagonal  and  five-diagonal  matrices. 

2.6.1  Sturm  sequence  for  tridiagonal  and  five-diagonal  matrices 

The  Sturm  sequence  is  a  method  for  calculating  the  number  of  roots  that  a  polynomial  has  in 
a  given  interval.  In  this  dissertation,  the  Sturm  sequence  method  for  band  matrices  is  used 
to  calculate  the  number  of  negative  eigenvalues  of  a  matrix.  The  following  theorems  can  be 
found,  for  example,  in  m  and  m- 
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Theorem  2.16.  Sturm  sequence  for  tridiagonal  matrices.  Let  A  be  a  symmetric  N  x  N 
tridiagonal  matrix,  and  let  Ak}k  =  a  1  where  k  =  1..IV,  Akjk+ 1  =  Ak+ijk  =  bk+\  where  k  = 
1..N  —  1.  All  other  elements  of  A  are  0. 

We  define  the  sequences  mk  and  qk  as 

m0  =1 

m\  =a\  (2.53) 

mk  =aimk-i  -  b2kmk-2  ,  k  =  2,3,...,N 

The  number  of  sign  changes  in  the  sequence  mk  is  the  number  of  eigenvalues  of  A  that  are 
smaller  than  0. 

Theorem  2.17.  Sturm  sequence  for  symmetric  five-diagonal  matrices.  Let  A  be  a 

symmetric  N  x  N  five-diagonal  matrix,  and  let  Ak>k  =  a\  where  k  =  1..N,  Ak>k+i  =  Ak+itk  = 
bk+\  where  k  =  1..N  —  1  and  Ak)k+ 2  =  Ak+2,k  =  ck+ 2  where  k  =  1..N  —  2. 

We  define  the  sequences  rnk  and  qk  as 

qk=  0  ,  k  <0 
mk  =0  ,  k  <  0 
m0  =1 

(2.54) 

ft_2  =6fc-imfc_3  -  Ci_i®_3  ,  k  =  3,4,  ...,N 
mk  =akmk-\  —  bkmk- 2  —  cl(ak-imks  —  ck_imk- 4)  +  2bkckqk-2  , 
k  =  1,2,..., AT 

The  number  of  sign  changes  in  the  sequence  mk  is  the  number  of  eigenvalues  of  A  that  are 
smaller  than  0. 

Remark  2.18.  In  implementations  of  this  method,  some  scaling  of  the  sequence  is  sometimes 
required  in  order  to  avoid  overflows  and  underflows  (see,  for  example,  pm)- 
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Suppose  that  we  wish  to  calculate  An,  the  n-th  largest  eigenvalue  of  the  tridiagonal  or 
five-diagonal  matrix  A.  Let  5  >  0.  We  observe  that  the  n-th  largest  eigenvalue  of  the  matrix 
(A  —  (An  +  5)1 )  is  negative.  Therefore,  the  number  of  sign  changes  in  the  sequence  m*,  for 
the  matrix  (A  —  (An  +  5)1)  is  no  smaller  than  n.  Similarly,  the  n-th  largest  eigenvalue  of  the 
matrix  (A  —  ( An  —  5)1)  is  positive.  Therefore,  the  number  of  sign  changes  in  the  sequence  m k 
for  the  matrix  {A  —  (An  —  5)1)  is  strictly  smaller  than  n. 

We  set  a  search  range  (aq,  02);  we  use  the  Sturm  sequence  to  verify  that  An  is  in  the  range, 
otherwise  we  extend  the  search  range.  We  then  use  bisection  to  narrow  the  range  (aq,  a^)  until 
«2  —  aq  is  smaller  than  the  desired  precision.  An  is  contained  within  the  range,  so  ( oq  +  a2)/2 
is  a  sufficient  approximation  for  An. 

2.6.2  The  inverse  power  method  for  tridiagonal  and  five-diagonal  matrices 

Let  B  be  a  symmetric  matrix,  and  let  An  7^  0  be  eigenvalue  of  B  with  the  largest  magnitude. 
Suppose  that  there  is  some  5  >  0  such  that  for  any  other  eigenvalue  Am  of  B ,  we  have  |An|  > 
(1  +  <5)|Am|.  The  power  method  is  a  well-known  method  for  calculating  the  eigenvector  v  and 
the  eigenvalue  An  by  iterative  calculation  of 

v(k+l)  =  Bv{k ).  (2.55) 

After  a  sufficient  number  of  iterations, 

Bv (fc)  «  Anv{k).  (2.56) 

Let  A  be  a  symmetric  tridiagonal  or  five-diagonal  matrix,  and  let  An  /  0  be  an  eigenvalue 
of  A  with  multiplicity  one.  Then  there  exists  5  >  0  such  that  for  any  other  eigenvalue  Am 
of  A,  |A„|(1  +  5)  <  |Am|.  The  inverse  power  method  is  a  well-known  method  for  calculating 
the  eigenvector  v  and  the  eigenvalue  An,  using  the  power  method  on  B  =  A-1.  Instead  of 
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computing  B  =  A  1  explicitly,  the  power  iteration  =  Bv^  is  computed  by  solving 


Av{k+1)  =w(fc). 


(2.57) 


2.6.3  Calculating  an  eigenvector  and  an  eigenvalue  of  a  tridiagonal  or  five- 
diagonal  matrix 

Let  A  be  a  symmetric  tridiagonal  or  five-diagonal  matrix.  Suppose  that  we  would  like  to 
calculate  the  n-th  largest  eigenvalue  Xn  and  the  corresponding  eigenvector  v  of  A,  such  that 

Av  =  Xnv  (2.58) 


Assume  that  Xn  has  multiplicity  one. 

First,  we  approximate  the  n-th  eigenvalue  Xn  using  the  Sturm  sequence  method  described 
and 


in  theorems 


2.16 


2.17 


We  require  the  approximation  Xn  to  be  close  to  Xn  compared  to  the 


difference  between  Xn  and  any  other  eigenvalue  of  A,  but  not  equal  to  Xn.  In  other  words: 


|  An  -  An  |  /  0 


(2.59) 


and 


|  An  An  |  <  |  An  Ani|  ,  Vm  ^  Tl 


(2.60) 


Next,  we  consider  the  matrix  ( A  —  X nI).  We  observe  that  the  eigenvector  v  that  we  wish 
to  calculate  is  also  an  eigenvector  of  ( A  —  X nI),  with  the  eigenvalue  an  =  Xn  —  Xn  /  0.  We 
observe  that  an  is  smaller  in  magnitude  than  any  other  eigenvalue  am  of  ( A  —  A nI).  We  use 
the  inverse  power  method  to  calculate  v. 


Finally,  we  obtain  a  better  estimate  for  eigenvalue  Xn  using  (2.58). 
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2.7  The  truncated  Laplace  transform 


Definition  2.19.  For  given  0  <  a  <  b  <  oo,  the  truncated  Laplace  transform  Ca $  is  a  linear 
mapping  L2(a,  b )  — >  L2( 0,  oo),  defined  by  the  formula 

(£a,b(f))  M  =  [  e-ta7(t)dt,  (2.61) 

J  a 

where  0  <  ui  <  oo. 

The  adjoint  operator  of  Ca ^  is  denoted  by  .  Obviously: 

COO 

((£a,b)*(g))(t)=  e~twg(u)  dw.  (2.62) 

■Jo 

The  operators  Ca,b  and  (Ca,b)*  are  compact  and  injective,  the  range  of  (£a,b)*  is  dense  in 
L2(o,  b )  and  the  range  of  Ca^  is  dense  in  L2( 0,  oo)  (see,  for  example,  ID- 


2.8  The  SVD  of  the  truncated  Laplace  transform 


In  this  section,  we  present  the  SVD  of  the  truncated  Laplace  transform,  which  is  the  main  tool 
we  use  to  investigate  the  properties  of  this  operator  in  this  dissertation. 

The  kernel  K  :  (0,  oo)  x  (a,  b)  — >  M  of  the  integral  operator  Ca ^  (defined  in  (2.61 ))  is  defined 
by  the  formula 


K(u,t)  =  e~wt, 


(2.63) 


so  that 


MS))  M  =  [b 

J  a 


K(u,t)f{t)dt. 


(2.64) 
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By  theorem  2.14,  there  exist  two  orthonormal  sequences  of  functions  un  E  L2(a,b )  and  vn  E 
L2(0,oo)  such  that 


K(u,t )  =  y ^vn(uj)snUn(t), 


(2.65) 


n= 0 


■^a,b(^n)  —  (^nvni 


(2.66) 


and 


(£-a,b)  (^n)  —  OinUn. 


(2.67) 


We  refer  to  the  functions  un(f)  as  the  right  singular  functions,  and  to  the  functions  vn(uj) 
as  the  left  singular  functions.  We  refer  to  an  >  0  as  the  singular  values.  The  functions  are 
numbered  n  =  0,1, ..,  and  they  are  sorted  according  to  the  singular  values,  in  descending  order. 

Observation  2.20.  The  multiplicity  of  an  in  this  decomposition  of  Ca,b  is  one  (see  0)- 

Observation  2.21.  A  simple  calculation  shows  that  (£a,b)*  o  £ab  of  a  function  /  E  L2(a,  b )  is 
given  by  the  formula 

(((£a,b)*  °£a,b)  (/))  (t)  =  [  r\—f(s)ds-  (2-68) 

J  a  t  i  S 


Clearly,  (£a,b)*  °  £a,b  is  a  symmetric  positive  semidefinite  compact  operator.  By  theorem 


2.14,  the  right  singular  functions  un  of  the  operator  Ca  b  are  also  the  eigenfunctions  of  the 


operator  (£a,&)*  °  £a,b,  and  the  singular  values  an  are  the  square  roots  of  the  eigenvalues  of 
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(. Ca,b )*  °  £ a,b ■  In  other  words, 


(((£a,b)*  °  £a,b)  ( Un ))  (t) 


t  +  S 


un(s)ds  =  a2nun{t). 


(2.69) 


Observation  2.22.  Similarly,  C,aj,°  ( [C-a,b )*  of  a  function  g  e  L2(0,oo)  is  given  by  the  formula 


poo  —a(ui+p)  ,  -b(u+p) 

ma,„  o  (A^n  ( 9 ))  m  =  /  - - p - »(p)dp. 

Jo  w  +  p 


(2.70) 


By  theorem  2.14,  the  left  singular  functions  vn  of  Ca ^  are  the  eigenfunctions  of  Ca,b  °  (£ a,b)* 
and  the  singular  values  an  are  the  square  roots  of  the  eigenvalues  of  £a,6  °  (£a,b)*-  In  other 
words, 


((£a,b  °  (£a,b)*)  (vn))  (w)  =  f 

Jo 


oo  p-a(u+p)  _|_  e-b(u+p) 


OJ  +  p 


-vn(p)dp  =  a2nvn(u).  (2.71) 


2.9  A  differential  operator  related  to  the  right  singular  func¬ 
tions  Ur 


Ln 


It  has  been  observed  in  [5J  that  the  integral  operator  (£a,b)*  °  (defined  in  (2.68))  commutes 
with  a  differential  operator. 


Theorem  2.23.  The  differential  operator  Dt,  defined  by  the  formula 


Dt(f)  )(t)  =  ft  (  (t2  -  a2)(b2  -  t2)p(f)  )  -  2{t2  -  a2) f{t), 


(2.72) 


commutes  with  the  integral  operator  (£a,b)*  °  £a,b  (defined  in  (2.68))  in  L2{a,b) 


It  has  also  been  shown  in  [5]  that  the  eigenvalues  of  the  operators  (Ca<b)*  °  £a,b  and  Dt  have 


a  multiplicity  of  one.  It  follows  from  theorem  2.23,  and  the  multiplicity  of  the  eigenvalues,  that 
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the  eigenfunctions  of  the  integral  operator  (Ca^)*  o  Ca  b  are  the  regular  eigenfunctions  of  the 


differential  operator  Dt.  By  (2.69),  these  functions  are  the  right  singular  functions  un  of  £a,b- 


Furthermore,  it  has  been  shown  that  if  the  eigenfunctions  of  Dt  are  sorted  according  to 
the  eigenvalues  of  Dt,  in  descending  order,  the  ra-th  eigenfunction  of  Dt  is  the  n-th  singular 
function  of  £a,fe-  Therefore,  un  is  both  the  n  +  1-th  right  singular  function  of  Ca,b,  the  n  +  1-th 
eigenfunction  of  (£a,b)*  o  £afe,  and  the  n  +  1-th  eigenfunction  of  Dt. 

We  denote  the  eigenvalues  of  the  differential  operator  Dt  by  xn.  By  theorem 
solution  to  the  differential  equation 


2.23 


un  is  the 


Dt(un ) )  (t)  =  ^  f(t2  -  a2)(b2  -  t2)^un(t)  )  -  2(f2  -  a2)un(t)  =  Xn,un(t). 


(2.73) 


2.10  The  function  i))n  associated  with  the  right  singular  function 

Un. 


The  right  singular  functions  un  of  £a  b  (the  operator  defined  in  (|2.61 ))  are  defined  on  the 
interval  (a,  b).  It  is  convenient  to  scale  and  shift  the  interval  (a,  b)  to  (0, 1). 

We  define  the  variable  x  G  (0, 1)  by  the  formula 


x  = 


t  —  a 
b  —  a 


t  =  a  +  (6  —  a)x. 


(2.74) 


The  functions  ipk  are  defined  using  the  change  of  variables  (2.74),  as  follows 


Definition  2.24.  The  function  ipn(x)  is  defined  via  the  corresponding  right  singular  function 
un,  by  the  formula 


=  Vb  —  a  un(a  +  (b  —  a)x). 


(2.75) 
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Observation  2.25.  Since  the  function  un  is  normalized  on  ( a,b ),  it  is  clear  from  (2.75)  that 
ipn  is  normalized  on  (0, 1) 


[  O n(x)Y 

■JO 


dx  =  1, 


(2.76) 


and  that  the  sequence  of  functions  ipn  forms  an  orthonormal  basis  in  L2( 0, 1). 


By  (2.68)  and  (2.74),  the  functions  ip n  are  the  eigenfunctions  of  the  integral  operator  T*  oT, 


where  T*  o  T  of  a  function  /  is  defined  by  the  formula 


T*oT)f)(x)  = 


1 


/  o  x  +  y  +  /3 


/(y)dy, 


(2.77) 


and  where  f3  is  defined  by  the  the  formula: 


P  = 


2  a 

b  —  a 


(2.78) 


Clearly,  T*  o  T  has  the  same  eigenvalues  as  ( Ca ^)*  o  Ca  b: 


((T*  o  T)  ^n)){x)  = 


1 


x  +  y  +  (3 


ipn(y)dy  =  a2n*pn{x) 


(2.79) 


Similarly,  by  (2.72)  and  (2.74),  ipn  are  the  eigenfunctions  of  the  differential  operator  Dx , 


which  is  defined  by  the  formula 


( Dx(f ))  (x)  =  fx(l  -  x)(/3  +  x)(/3  +  1  +  x)^/(x)  )  — 2x(x+/3)/(x).  (2.80) 
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In  other  words 


(Dx(ipn))  (x)  = 

(x(l  -  x)(/3  +  x)(/3  +  1  +  x)-^^n(x)^j  -  2x(x  +  p)'ipn(x)  =  (2.81) 

Xnlpn(x), 


where  Xn  are  the  eigenvalues  of  Dx. 


2.11  A  differential  operator  related  to  the  left  singular  func¬ 
tions  Vn 


It  has  been  observed  in  [5j  that  the  integral  operator  Ca o  (Ca,b)*  (defined  in  (2.70))  commutes 
with  a  differential  operator. 

Theorem  2.26.  The  differential  operator  Dw,  defined  by  the  formula 


DM))  M  =  [[Ca,boDo(jra>b) 
d2  ' 


-l 


M  = 


dw2 


+  (a2  +  62)^  ^w2-j^/(w)^  +  (- a2b2u 2  +  2a2)  f(u), 


(2.82) 


commutes  with  the  integral  operator  £a,6  °  (£a,b)*  ( defined  in  (2.70)).  The  left  singular  func¬ 
tions  vn  are  the  eigenfunctions  of  Du . 


We  denote  the  eigenvalues  of  D u  by  x\-  By  theorem 


2.26 


the  function  vn  is  the  solution 
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of  the  differential  equation 


(A^(ffc))  M  = 

=  (^Tj2 yfc(w))  +  («2  +  +  (_a2&2cj2  +  2«2)Dc(w) 


=  Xfc^fc(w). 


(2.83) 


Observation  2.27.  The  eigenvalues  of  are  equal  to  the  eigenvalues  of  Dt: 


Xn  =  X*n 


(2.84) 


2.12  The  functions 


Having  introduced  the  operator  Ca j,  in  (2.61)  and  its  adjoint  (£a,b)*  in  (2.62),  we  now  discuss 


the  properties  of  the  function  generated  by  applying  (£a,b)*  to  the  Laguerre  function  (defined 


in  (2.47)).  By  (2.42),  (2.47)  and  (2.61) 


coo  COO 

((£„,&)*  ($*))  (t)  =  /  e~ut$k(u,)du,  =  /  e~^t+1^Lk(u)du  = 

Jo  Jo 


0 

-k- 1 


(2.85) 


In  particular,  at  t  =  1/2,  (2.85)  becomes 


COO 

(M*ft))(l/2)=  /  e-«/2<l>fc(g)dg 

4o 


= 


1  if  k  =  (1 
0  otherwise 


(2.86) 


Differentiating  (2.85),  we  obtain 


((Cat)*  ($*))' (t)  =  (8 k  -8 1  +  4)(2t  -  l)fe-1(2t  +  l)-fc-2,  (2.87) 
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which,  at  t  =  1/2,  becomes 


((£«,br  (**))' (1/2) 


—  1  if  k  =  0 

\  1  if  k  =  1 


0  otherwise 


(2.88) 
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Chapter  3 


Analytical  apparatus 


In  this  part  of  the  dissertation,  we  discuss  certain  useful  properties  of  the  truncated  Laplace 
transform;  we  begin  with  a  brief  discussion  of  the  scaling  properties  of  the  truncated  Laplace 
transform,  and  with  a  definition  of  a  standard  form  of  the  truncated  Laplace  transform.  We 
proceed  to  define  the  transform  C7  and  discuss  various  symmetry  properties  associated  with 
it.  We  then  discuss  the  expansions  of  un  and  vn  in  orthonormal  bases,  and  show  that  the 
calculations  of  un  and  vn  can  be  phrased  as  benign  eigensystem  calculations.  This  chapter  is 
concluded  with  brief  discussions  of  several  miscellaneous  useful  properties. 


3.1  On  the  scaling  properties  of  the  truncated  Laplace  trans¬ 
form 


The  truncated  Laplace  transform  (as  defined  in  ( 2.61[) )  can  be  generalized  to  the  form 


(3.1) 


with  arbitrary  0  <  c  <  oo,  0  <  a  <  b  <  oo. 
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Observation  3.1.  The  properties  of  the  truncated  Laplace  transform  are  determined  by  the 
ratio 


7  =  b/a  >  1 


(3.2) 


(see,  for  example,  [I]). 


Observation  3.2.  The  particular  choice 


a  =- 


2^ 

6=^ 

2  ’ 


c  =1, 


(3.3) 


yields  several  useful  properties,  which  we  will  discuss  in  this  dissertation. 


Due  to  observations  3.2  and  3.1,  in  the  remainder  of  this  dissertation  we  will  be  assuming 


without  loss  of  generality  that  the  values  of  a,b  and  c  are  as  defined  in  (3.3).  In  other  words, 


we  will  restrict  our  attention  to  the  following  form  of  the  truncated  Laplace  transform: 

Definition  3.3.  For  a  given  1  <  7  <  00,  we  will  denote  by  £7  :  yy)  — >  L2(0,oo)  the 

operator  defined  by 


Fry  -  F  1  77  -  F  y 


■rry  -  i—  1  77 

277’  2 


v7 
277’  2  ’ 


(3.4) 


The  operator  will  be  referred  to  as  the  ‘‘standard  form”  of  the  truncated  Laplace  transform. 
Obviously,  £7  of  a  function  /  €  L2(  ^=,  yy)  is  defined  by  the  formula 


(Fry(f))  (w)  =  77  2(/))(w)  =  /  2  e~tuf(t)dt. 

2v7>  2  d  J_^ 


(3.5) 
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Where  there  is  no  danger  of  confusion,  we  write  C  instead  of  £7,£a  ;,  and  £a,b,e>  and  we 


denote  the  adjoint  of  C  by  £*. 


Remark  3.4.  Combining  (3.2)  with  (2.78),  we  observe  that  the  quantity  (3  (defined  in  (2.78)) 
is  related  to  7  by  the  formula 


2  a 

b  —  a 


(3.6) 


Remark  3.5.  Let  un,  vn  and  an  be  the  n  +  1-th  right  singular  function,  left  singular  function 
and  singular  value  of  £~  ^  such  that 


—  OL'nVn-  (3.7) 

Let  7  =  b/ a  and  let  un,  vn  and  an  be  the  n  +  1-th  singular  functions  and  singular  value  of  £7. 
Then,  the  SVD  of  C~  r  ~  is  related  to  the  SVD  of  the  standard  form  £7  by: 

(■ ta/a ),  (3-8) 


UnM 


vn(uca/a), 


(3.9) 


and 


OLn 


—  &n/ \/c* 


(3.10) 
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3.2  The  transform  C1 

In  this  section  we  define  the  transform  C7  which  is  useful  in  the  discussion  of  certain  symmetry 
properties. 

Definition  3.6.  We  define  the  new  variable  seR  via 


s  =  21og(2f)/log(7). 


(3.11) 


For  7  >  1,  we  define  the  transform  C1  of  a  function  f  by  the  the  formula 


(C7(/))(»)=7‘/4/(t*/72). 


(3.12) 


Observation  3.7.  A  simple  calculation  shows  that 


fs(*2)  4  rt2 

/  (C^yU))  (s)  (C^(g))  (s)ds  =  — —  /  f(t)g(t)dt. 

Js(ti)  log  7  Jt,. 


(3.13) 


We  are  particularly  interested  in  the  case  where  a,  b  are  as  defined  in  (3.3).  In  this  case,  C7 
becomes  a  mapping  L 2  (2(77’  ^f)  (—1,  !)• 


3.2.1  The  functions  {C1  o  C*)  ($&) 


In  this  subsection,  we  discuss  certain  properties  of  the  Laguerre  functions  (defined  in  (2.47)), 
related  to  the  operator  C1  (defined  in  ( 3. 12[) ) . 


A  simple  calculation  shows  that  C7  of  the  function  £*(4>fc)  (see  (2.85))  is  given  by 


((C7  o  C*)  (<&*))  (s)  =  7S/4  (7S/2  -  l)"  (7S/2  +  l) 


-fe-i 


(3.14) 
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Clearly, 


|((C7o  £*)(<£,))  (s)|  = 


7 


s/ 4 


js/2  _|_  ^ 


7' 


s/2  _  l 


7 


s/2  _|_  l 


(3.15) 


from  which  it  immediately  follows  that 


|((Cyo£*)  (**))(s)|< 


1 


yS/2  _  ! 


s/2  +  1 


(3.16) 


Observation  3.8.  For  all  1  <  7  <  00  and  s  6  R,  we  have 


s/2! 


s/2  +  l 


<  1;  it  is  therefore  obvious 


from  (3.15)  that  |  ((C7  o  £*)  (3>j.))  (s)|  decays  exponentially  as  k  grows. 


Observation  3.9.  By  (3.14), 


((C7  o  £*)  (*fc))  (a)  =  (-l)fc  ((C7  o  C* )  ($*))  (-a). 


(3.17) 


In  other  words,  for  an  even  k,  the  function  (C7  o  £*(<&&))  ( s )  is  even;  and  for  an  odd  k,  it  is 
odd. 


Observation  3.10.  By  (3.14),  at  the  point  s  =  0, 


((cl  o  r)  (<&*))  (o)  = 


1/2  if  fc  =  0 
0  otherwise 


(3.18) 


Observation  3.11.  By  differentiating  (3.14)  and  setting  s  =  0  we  obtain: 


((Gy  o  £*)($*))'  (0)  =  < 


l°g(7)/4  if  k  =  1 
0  otherwise 


(3.19) 
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3.2.2  C7  of  the  right  singular  function  un 

Definition  3.12.  We  introduce  the  function  Un,  which  we  define  by  the  formula 


Un(s)  =  (C7(un))  (s), 


(3.20) 


where  un  is  a  right  singular  function  of  the  operator  Ca,b,  and  C7  is  defined  in  (3.12). 


By  (3.13),  and  since  un  is  normalized  on  (a,  b),  the  norm  of  Un  on  ^2 ,  2 ^  is: 


2log(2 b) 

c  iog(7)  q  4 

/  (Un(s))2ds  =  - - 

/2log(2a)  V  logy 

log(7) 


(3.21) 


Equation  (3.21 )  holds  for  an  arbitrary  choice  of  a  and  b  such  that  b/a  =  7.  In  this  dissertation 
we  assume  a  =  7^7,6  =  77  (as  defined  in  (3.3)).  By  substituting  (3.3)  into  (3.21),  the  interval 
s  G  (2To@rP  2lSf )  becomes  0  G  (-1, 1). 


In  the  case  of  £7  (the  standard  form  of  the  truncated  Laplace  transform,  defined  in  (3.5)), 


by  (2.72)  and  (3.20),  the  functions  Un  are  the  eigenfunctions  of  the  differential  operator  Ds, 


defined  by  the  formula 


D.«(f) )  (s)  =  (log  (^7))  2  ^  (72  +  1  -  2y  cosh  (2s  log  (-/y)))  ^/(s) 


ds" 


1 


7 


-  (  2^ cosh(2'sl°g(v/7))  +  47  -  4  )/(«)■ 


(3.22) 


A  simple  calculation  shows  that  the  eigenvalues  of  Ds,  which  we  denote  by  p,n,  are  related  to 


the  eigenvalues  Xn  (defined  in  (2.72))  by  the  formula: 


hn  =  47Xn. 


(3.23) 
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3.3  The  symmetry  property  of  un  and  JJn 


By  0,  the  right  singular  functions  un  of  Laf,  (the  operator  defined  in  (|2.61[))  satisfy  a  form 


of  symmetry  around  the  point  \fab.  In  the  case  of  the  standard  form  £7,  defined  in  (3.5),  we 
have  \fab  =1/2,  and  the  symmetry  relation  is: 


un  (  ^  )  =  (— 1)”2 tun(t) 


(3.24) 


Observation  3.13.  In  the  case  of  standard  form  £7,  it  follows  from  (3.20),  that  the  functions 


Un  (defined  in  ( 3.20|) )  are  even  and  odd  functions  in  the  regular  sense: 


Un(s)  =  (C7(un))  (S)  =  (-1  )nUn(s). 


(3.25) 


In  particular,  at  the  point  s  =  0,  we  have: 


U2j+l(0)  —  (Cj(u2j+ 1))  (0)  —  0, 


(3.26) 


and 


UhjiO)  =  (Cy(u2j))'  (0)  =  0. 


(3.27) 


Remark  3.14.  The  functions  Un  are  even  and  odd  functions  around  the  point  s  =  0  in  the 


case  of  the  standard  form  £7  (as  defined  in  (3.5)).  Similar  symmetry  exists  for  C7  of  the 


right  singular  functions  of  Ca ^  (as  defined  in  (2.61[)),  however  the  center  of  symmetry  is  not 
necessarily  s  =  0. 
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3.4  The  differential  operator  Dx  and  the  expansion  of  t(jn  in  the 


basis  of 


In  this  section  we  consider  the  expansion  of  functions  /  E  L2(0, 1)  in  the  orthonormal  basis  of 


the  polynomials  Pf  (defined  in  (2.16)): 


/(*)  =  Y2hkPpx). 


(3.28) 


k= 0 


Lemma 


3.15 


describes  the  operation  of  Dx  (defined  in  (2.80))  on  a  basis  function  Pj*.  The 


result  is  used  to  express  the  functions  ifn  (dehned  in  (2.75))  via  a  five-terms  recurrence  relation 


or  a  solution  to  a  benign  eigensystem,  specified  in  theorem  3.16 


Lemma  3.15.  Applying  the  differential  operator  Dx  to  the  polynomial  Pf  yields  a  linear  com¬ 
bination  of  Pf_2,P*_vPf,P;+}  and  P*+ 


( Dx(P*k ))  (*)  = 

— _ (■ k—l)2k ~ 

4^2fc^3(2fc-l)x/2fc+T  k- 

_  fc3(l+/3)  p*  /  \ 

v/2fc4Tv/2fc+T  k-iy X) 

(— 4—6/3— 2fc/3(2+3/3)+fc2(7+12,3+2/32)  +  (2fc3+fc4)(7+16/3+8/32))  -p*/ 

2(2fc-l)(2fc+3)  Pk  \X 

(fc+l)3(l+/3)  D*  ( ^ 
v/2FPTv/2FP3  k+lV  ) 

(fc+l)2(fc+2)2  p* 

4  v/2fc+T(2fc+3)  ^2fc+5  k+2  V  >  > 


where  B  =  Tp- 

“  b—a 


^j-  (as  defined  in  3.6). 


Proof.  By  the  definition  of  Dx  (in  (2.80)), 


(Dx(P*k))(x)  = 

=  ^  ({P  +  x){pP  1  +  x)x(l  -  -2x(x  +  P)Pf(x). 


(3.30) 
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Using  the  chain  rule, 


(ArTO)  (*)  = 

=  (^(/3  +  X)(/?  +  1  +  X))  _X)^Pfc*(X)) 

+  (/3  +  x)(/3  + 1 +  x)-|^  (x(l  -  x)^jPfc(x))  -  2x(x  +  /3)Pfe*(x)  = 

=(1  +  2x  +  2/3)  (x(l-x)-^Pfe*(x)) 

+  (x2  +  x(l  +  2/3)  +  /3  +  /32)  —  ^x(l  —  x)— P ^  (x)^ 


-  2x(x  +  /3)P£(x) 


Using  identities  (2.12),  (2.13)  and  (2.14), 


(Dx(P*k))(x)  = 

(-i  +  fc)2  ep*_2(x) 

4(-l  +  2/s)(l  +  2k) 
fe3(l  +  ^)P*_i(x) 

1  +  2fc 

(— 4+7fc2+14fc3+7fc4— 6/3— 4fc/3+12fc2/3+32fc3/3+16fc4/3— 6fc/32+2fc2/32+16fc3/32+8fc4/3 

2(-l  +  2fe)(3  +  2fc) 

(l  +  fc)3(l  +  /3)P*+1(x) 

1  +  2fc 

(l  +  fc)2(2  +  fc)2P*+2(x) 

4(1 +  2fc)(3  +  2fc) 


Finally,  substituting  (2.11)  into  (3.32)  gives  (3.29) 


Theorem  3.16.  Fe£  t/ie  function  ipn(x)  be  as  defined  in  (2.75).  Let  hn  = 


(3.31) 


21 

APk{x) 

(3.32) 

□ 

.)T  be  the 
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vector  of  coefficients  in  the  expansion  of  i/jn(x)  in  the  basis  of  the  polynomials  Pjf: 


^ n{x )  =  ££ kPk(x) 
k= 0 


(3.33) 


Then,  hn  is  the  n  +  1-th  eigenvector  of  M  : 


Mhn  =  Xnhn , 


(3.34) 


where  M  is  the  five-diagonal  matrix 


Mk-2,k  = 


(fc-i)2fc2 

4v/2fc=3(2fc-l)v/2FFT 


Mk-l,k 

Mk,k 


fc3(l+/3) 

v/2fc^Tv/2fc+T 

(— 4—6/3— 2fc/3(2+3/3)+fc2  (7+12/3+2/32  )  +  (2fc3+fc4)  (7+16/3+8,32  )  ) 
2(2fc— l)(2fc+3) 


Mi 


k+l,k 


 (fc+l)3(l+/3) 


^/2FFT^/2FF3 


Mk+2,k  = 


(fc+l)2(fc+2)2 
4v/2fc+T(2fc+3)\/2FF5  ’ 


(3.35) 


and  where  \n  are  the  eigenvalues  of  the  differential  operator  Dx,  and  k  =  0, 1,  2.... 


Proof.  By  (2.81),  ipn(x)  is  an  eigenfunction  of  Dx,  with  the  eigenvalue  \n.  Since  the  differential 


operator  is  linear, 


(Actyn))  (x)  =  J2K  (. DX(P *))  (x)  =  Xn  £  KP*(x). 

k= 0  /c=0 


(3.36) 


Using  lemma  3.15,  (3.36)  becomes  (3.34). 


□ 


Observation  3.17.  Clearly,  hf  is  the  inner  product  of  ifn  and  Pf : 


hk=  Pf{x)i/jn(x)dx. 

Jo 


(3.37) 
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3.4.1  The  decay  of  the  coefficients  in  the  expansion  of  ipn 


Since  the  functions  'ipn  are  smooth  regular  solutions  of  a  differential  operator,  they  can  be 
efficiently  expressed  using  an  orthogonal  basis  of  polynomials.  In  other  words,  we  expect  the 
coefficients  in  the  expansion  of  in  terms  of  the  polynomials  Pk  to  decay  rapidly.  In  this 
subsection,  we  obtain  a  bound  for  this  decay. 

Lemma  3.18.  Let  0  <  (5  <  oo  and  0  <  y  <  1.  We  introduce  the  notation 

P  =  7(1  +  (2/5)2)  -  1  =  VW+T).  (3.38) 


Then, 


i 


1  ^ 


x  +  y  +  j3 


P£(x)dx 


2 

d  y  < 


( 

V 


2y/2k  +  1 


(3.39) 


where  P£  is  defined  in  (2.16). 


Proof.  We  recall  from  (2.24)  that 


(■ x  +  y  +  fi )  1P£(x)dx 


=  2Q*k{y  +  13  +  l)72^Tl, 


(3.40) 


where  Q*k  is  defined  in  (2.23).  So,  by  lemma  2.8 


(x  +  y  +  j3)  1P?(x)dx 


< 


2y/2k  +  1 
fc+i 

1  + 


log  2 


1+7 

P  . 


+  1 


(3.41) 


By  squaring  (3.41)  and  integrating  over  y,  we  obtain  (3.39) 


□ 


Lemma  3.19.  Let  hk  be  the  k  +  1-th  coefficient  in  the  expansion  defined  in  (3.33),  of  the 
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function  ifn  ( defined  in  (2.75))  in  the  basis  of  the  polynomials  Pf  ( defined  in  (2.16)).  Then, 


— 


2 2^2k  +  l 


-\  fc+i 

1  +  /?) 


log  2 


1  +  /3 


+i 


(3.42) 


Where 


0=^(1  +  (2/3)2)  _  1  =  v/4/3(l  +  /3) 


(3.43) 


and  (5  is  as  defined  in  (3.6). 


Proof.  We  substitute  (2.79)  into  (3.37)  and  change  the  order  of  integration 


K  =  an2 


1  r  1 


Pf(x)^n(y)dxdy  = 


=  a, 


o  Jo  x  +  y  +  P 
-([  «»>(/ 


(3.44) 


By  the  Cauchy-Schwarz  inequality, 


n|  <«"2i 


L 


2dyi/  /  (  f  -i  „  1  ,  )  dy- 


o  Vjo  x  +  y  +  P 


(3.45) 


By  (2.76)  and  (3.39) 


/ 


Ityfel  < 


h  2T  +1 


U1  +  ^ 


(3.46) 


□ 
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3.5  The  differential  operator  Du  and  the  expansion  of  vn  in  the 
basis  of 

In  this  section  we  consider  the  expansion  of  functions  g  6  L2( 0,  oo)  in  the  basis  of  the  Laguerre 
functions  (the  functions  defined  in  ( 2.47|) ) : 


sM  =  'Y^r}k$k(u). 


(3.47) 


k= 0 


Lemma  3.21  describes  the  operation  of  the  differential  operator  Du  (define  in  (2.83))  on 
This  relation  is  used  to  express  the  expansion  of  the  left  singular  function  vn  of  the  operator 


j£a,b  (the  operator  defined  in  (2.61))  via  a  five-terms  recurrence  relation,  or  as  a  solution  to  a 


benign  eigensystem  described  in  theorem  3.22 


Remark  3.20.  Lemma  3.21[  theorem  |3.22  and  the  discussion  in  section  3.5.5  apply  to  the 


operators  associated  with  Ca,b  (defined  in  (2.61))  with  an  arbitrary  choice  of  0  <  a  <  b  <  oo. 
Subsections 


3.5.1 


and 


3.5.2 


apply  to  the  special  cases  of  £  172,7/2  and  £1/(27), 1/2-  Subsections 


3.5.3  and  3.5.4  treat  to  the  standard  form  of  the  truncated  Laplace  transform  £~,  as  defined 


in  (3.5). 


Lemma  3.21.  Applying  the  differential  operator  Du  (defined  in  (2.83))  to  the  Laguerre  func 


tion  &k  (defined  in  (2.J^1))  yields  a  linear  combination  of  the  Laguerre  functions  &k-2,  &k-i, 
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®k,  $fc+i  and  $k+2: 

(£>U $*))  (w)  = 

-  ^  (4«2  -  !)  (4&2  -  !)  (&  -  l)fc^fc-2  (w) 

+  |/c2  (l6a262  -  l)  $fc_i  (w) 

+  ^  (fc(fc  +  1)  (—48 a2b2  —  4 a2  —  4b2  —  3)  +  (— 16a2 62  +  12a2  —  462  —  l))  3>*.  (oj) 
+  +  l)2  (l6a26“  —  l)  (a;) 

_  16  ^°2  ~~  (4^  ~~  2)(fc  +  l)<I,fc_|_2  (w) . 

(3.48) 


Proof.  Applying  Du  to  a  Laguerre  function  yields 

(pu{*kj)  (®)  = 

=  +  (°2  +  fe2)^w2^$fc(:c)  +  (- a2b2u] 2  +  2a2)$fc(w) 

(3.49) 


A  somewhat  tedious  derivation  from  (3.49),  using  identities  (2.39),  (2.40)  and  (2.41),  yields 

□ 


(3.48). 


Theorem  3.22.  Let  vn{ui)  be  the  n  +  1-th  left  singular  function  of  the  truncated  Laplace 
transform.  Let  rf1  =  (rfi ,  rjf,  ...)T  be  the  vector  of  coefficients  in  the  expansion  of  vn(u)  in  the 


basis  of  Laguerre  functions  (the  functions  defined  in  (2.4'))),  such  that 


vn(uj)  =  ^2ri^k(u). 


(3.50) 


k= o 
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Then,  rjn  is  the  n  +  1-th  eigenvector  of  M : 


MVn  =  xW 


(3.51) 


where  M  is  the  symmetric  five- diagonal  matrix 

Mk-2,k  =  -  ^  (4a2  -  1)  (4 b2  -  1)  (k  -  1  )k 
Mk-i.k  =  +  \k2  (l6a262  -  l) 

Mk  k  =  +  i  {k{k  +  1)  (— 48a262  -  4a2  -  462  -  3)  +  (-16aV  +  12a2  -  462  -  l)) 
8 

Mk+i,k  =  +  +  l)2  (l6a262  —  l) 

Mk+2,k  =  -  ^  (4a2  -  1)  (462  -  1)  (k  +  2 ){k  +  1), 

(3.52) 


and  where  %*  are  the  eigenvalues  of  ( defined  in  (2.83)),  and  k  =  0, 1,2.... 


Proof.  By  (2.83),  the  left  singular  function  vn  is  an  eigenfunction  of  the  differential  operator 


Du,  and  therefore 


A^n))  M  =  x*vn{u). 


(3.53) 


Substituting  (3.50)  into  (3.53)  and  using  the  linearity  of  the  differential  operator,  we  obtain: 


£  ni  (£U it)j  M  =  x"  £ 

k= 0  k= 0 


(3.54) 


Using  lemma  3.21  and  (3.54),  we  obtain  (3.51). 


□ 
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Observation  3.23.  Clearly, 


/•OO 

Vk=  vn(u)®k(u)du ;. 

Jo 


(3.55) 


Remark  3.24.  Expressing  the  left  singular  functions  vn  in  a  similar  way,  using  Hermite  polyno¬ 
mials  or  parabolic  cylinder  functions,  yields  a  similar  framework,  with  a  seven-diagonal  matrix 

M. 


3.5.1  A  special  case  of  theorem  3.22  a  =  1/2 


We  observe  that  there  are  two  special  choices  of  a  and  b  for  which  the  matrix  M  (defined  in 


(3.52))  becomes  tridiagonal.  We  briefly  describe  these  two  cases  in  this  subsection  and  in  the 
next  subsection. 


The  substitution  of  a  =  1/2,  b  =  7/2  into  (3.52)  yields  the  first  tridiagonal  case  of  M: 


Mk-i,k  =  +  ^  (V  -  l)  k2 

Mk,k  =  +  \  (-72  -  2  (72  +  1)  k2  -  2  (72  +  1)  k  +  1) 
Mk+l,k  =  +  "^  (72  —  l)  (&  +  l)2 


(3.56) 


3.5.2  A  special  case  of  theorem  3.22  6=1/2 


27’ 

Mk-l,k  = 
Mk,k  = 
Mk+l,k  = 


A  substituting  ofa  =  A,  b  =  1/2  into  (3.52)  yields  the  second  tridiagonal  case  of  M : 


(72  —  l)  k 2 
472 

(2  (72  +  l)  k2  +  2  (72  +  l)  k  +  72  -  l) 


472 


(3.57) 


(72  -  1)  (fc  +  if 
472 
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3.5.3  A  special  case  of  theorem  3.22  the  standard  form  of  the  truncated 


Laplace  transform,  as  defined  in  3.5 


We  now  consider  theorem  3.22  in  the  case  of  the  standard  form  £7  (defined  in  (|3.5|));  in  other 
words,  we  set  a  =  and  b  =  ^  (as  defined  in  (2.61)).  We  will  show  that  in  this  case,  the 


even-numbered  left  singular  functions  V2j  are  expanded  using  only  the  even-numbered  Laguerre 
functions  <3?2m,  and  that  the  odd-numbered  left  singular  functions  V2j+i  are  expanded  using  only 
the  odd-numbered  Laguerre  functions  $2m+i-  Furthermore,  we  will  show  that  the  expansions 
of  V2j  and  V2j+i  can  be  obtained  from  two  benign  tridiagonal  eigensystems. 


Observation  3.25.  We  substitute  a  =  j^,b  =  ^  (as  specified  in  (3.3))  into  (3.50).  We 
observe  that  the  first  off  diagonal  of  M  vanishes,  but  the  second  off  diagonal  does  not  vanish: 


Mk-2,k  ~  + 


Mk,k  ~  + 


Mk+2,k  =  + 


(7  —  1  )2{k  —  1  )k 
167 

((~72  -  67  -  l)  k(k  +  1)  -  72  -  27  +  3) 
87 

(7-  l)2(fc  +  l)(fc  +  2) 

I67 


(3.58) 


Let  Mj  j  be  an  entry  of  M  that  does  not  vanish.  Then,  we  observe  that  both  i  and  j  must 
be  even  or  both  must  be  odd.  In  other  words,  the  non-zero  elements  can  be  found  only  in  even- 
numbered  columns  of  even-numbered  rows,  and  in  odd-numbered  columns  of  odd-numbered 
rows  of  M. 

We  split  the  matrix  M  into  two  matrices;  one  of  the  matrices  contains  all  the  even  rows  of 
all  the  even  columns,  and  the  other  matrix  contains  all  the  odd  rows  of  all  the  odd  columns. 
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These  are  the  two  tridiagonal  matrices  Meven  and  Meven ,  specified  by  the  formulas: 


and 


71  /reven 


7i /reven  _ 

'  lj-j  ~ 


7+1,7 


(7-1)2(2j-1)2j 

I67 

((— 72  ~  67  -  1)  2j(2j  +  1)  -  72  -  27  +  3) 
87 

(7  —  l)2(2j  +  l)(2j  +  2) 

I67 


m: 


odd  _(7-l)2(2j)(2j  +  l) 


7-1,7 


I67 


Modd  _  ((— 72  —  67  —  l)  (2 j  +  l)(2j  +  2)  -  72  -  27  +  3) 


m: 


7,7  87 

(7  —  l)2(2j  +  2)  (2j  +  3) 


odd 


7+1,7 


167 


(3.59) 


(3.60) 


We  introduce  the  notation  7jeven,J  and  y*’et’e,t  for  the  j  +  1-th  eigenvector  and  eigenvalue  of 
Meven ,  and  rjodd  J  and  x*j’°dd  for  the  j  +  1-th  eigenvector  and  eigenvalue  of  Modd- 


j^even^evenj  _  ,even ^even,j 


=  Xi 


V 


(3.61) 


and 


M°ddiqOdd,j  _  ^*,odd^odd,j 


(3.62) 


Observation  3.26.  Let  y*  be  an  eigenvalue  of  M.  Then,  x*  is  either  an  eigenvalue  of  Meven , 
or  an  eigenvalue  of  Modd.  Any  eigenvalue  of  Meven  or  Modd  is  an  eigenvalue  of  M. 

Observation  3.27.  The  vector  (r)QVen^ ,  0,  r)f  ’en^ ,  0,  ....^  is  an  eigenvector  of  M  with  the 

1  *,even 

eigenvalue  \j 


Observation  3.28.  The  vector 


T 


is  an  eigenvector  of  M  with  the 
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i  *,odd 

eigenvalue  %  ■’ 

We  introduce  the  notation  VjVen(u ;),  v°dd(uj ): 

OO 

^enH  =  Er?renJ'^M  (3-63) 

1=0 


V°dd(uj)  =J2rfdd,:$2l+l{u). 
1=0 


(3.64) 


Observation  3.29.  Each  function  in  the  sequences  Vjven  and  v°dd  is  a  left  singular  function. 
Each  left  singular  function  is  either  in  the  sequence  of  functions  VjVen,  or  in  the  sequence  vodd . 


It  remains  to  be  shown  which  function,  in  which  of  the  two  sequences  VjVen  and  vodd , 
corresponds  to  the  n  +  1-th  left  singular  function  vn. 

Lemma  3.30.  Let  jfven^  be  the  j  +  1-th  eigenvector  of  Meven  and  let  v)ven  be  as  defined  in 


(3.63). 


Let  rfdd']  be  the  j  +  1-th  eigenvector  of  Modd  (defined  in  (3.60))  and  let  v)dd  be  as  defined 


in  (3.6f). 
Then, 


v2j(u)  = 

v2j+i{uj)  =  v°dd{u). 


(3.65) 


Proof.  By  observation  3.29,  v(jven  is  a  left  singular  function.  Let  um  be  the  corresponding  right 


singular  function  of  the  truncated  Laplace  transform  C  (the  operator  defined  in  (3.5)).  By 


(3.63)  and  (2.67), 


Um.  — 


_  „-i  ,  even,j 


m  /  v 
1=0 


even.)  /  r  \ 

Vi  (£7)  ($2 1)- 


(3.66) 
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We  multiply  both  sides  of  (3.66)  by  the  operator  C7  (as  defined  in  (3.12)),  and  use  the  definition 
of  Un  in  (3.20),  to  obtain 


Um  =  a'1  r,rnJ  (C7  o  (£7)* )  ($2{). 


(3.67) 


1=0 


By  (3.17),  the  functions  ((C7  o  (£7)*)  (<&2z))  (s)  are  even  functions,  so  Um  is  an  even  func¬ 


tion;  therefore,  by  (3.25),  m  must  be  an  even  number.  In  other  words,  v(ven  is  the  even- 


numbered  left  singular  function  vm. 

By  a  similar  argument,  v jdd  is  an  odd-numbered  left  singular  function.  In  other  words,  the 
sequence  of  functions  v?ven  is  the  sequence  of  even-numbered  left  singular  functions  vn  and  the 
sequence  of  functions  v°dd  is  the  sequence  of  odd-numbered  left  singular  functions  vn.  Based 


on  these  facts,  it  is  a  matter  of  simple  bookkeeping  to  obtain  (3.65)  using  observation  3.26 


□ 


3.5.4  Additional  properties  of  vn  in  the  case  of  the  standard  form  of  the 
truncated  Laplace  transform 


t)q  and  77”,  the  first  two  coefficients  in  the  expansion  (3.50)  of  vn,  are  related  to  Un{ 0)  (the 


function  defined  in  (3.20),  at  the  s  =  0)  and  to  the  value  of  the  derivative  U'n{ 0). 


Lemma  3.31.  Let  un  and  vn  be  the  n  +  1-th  right  and  left  singular  function  of  £7  (defined  in 
(3.5)).  Let  rjO  and  rff  be  the  first  and  second  coefficients  in  the  expansion  defined  in  \3.5Cfy,  of 


vn  in  the  basis  of  Laguerre  functions  (the  functions  defined  in  (2.  f.1)).  Let  Un  be  C7un,  as 


defined  in  (3.20).  Then: 


Un{  0)  =  a~lr,Z/ 2, 


(3.68) 
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and 


K(  0)  =  olu  V  log(7)/4. 


(3.69) 


Proof.  By  (3.50)  and  (2.67), 


un  = 


k= 0 


(3.70) 


We  apply  the  operator  C7  (defined  in  (3.12))  to  (3.70),  and  use  (3.20)  to  obtain 


Un(s)  =  a'1  ((C7  o  (£7)*)  ($,))  (S) 


k= 0 


In  particular,  at  the  point  s  =  0, 


14(0)  =  a”1  ]T  Vk  ((C't  °  (£,)*)  ($fc))  (0) 


fc=0 


(3.71) 


(3.72) 


We  then  use  (3.18)  to  obtain  (3.68) 


We  differentiate  (3.71)  and  set  s  =  0 


Km  =  K1  E  Vi  ((A ° (A)*)  (*U)'  (0). 


fc=0 


(3.73) 


We  use  (3.19)  to  obtain  (3.69). 


a 


Remark  3.32.  Similar  relations  for  the  value  of  the  right  singular  function  itn(l/2)  of  £7  (3.5) 


at  t  =  1/2  and  for  the  derivative  u'n(  1/2)  are  easy  to  obtain  from  lemma  3.31  or  by  a  similar 
construction. 


Remark  3.33.  In  the  other  spacial  cases  of  Ca,b >  where  a  =  1/2  or  b  =  1/2,  similar  rela¬ 
tions  exist  between  the  value  of  the  function  un  at  the  ends  of  the  interval  [a,  b]  and  the  first 
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coefficients  in  the  expansion. 


3.5.5  Decay  of  the  coefficients  in  the  expansion  of  vn  in  the  basis  of  <E>fc 
The  left  singular  functions  vn  are  smooth  functions,  and  they  are  therefore  efficiently  expressed 


using  Laguerre  functions  (the  functions  defined  in  (2.47)).  In  other  words,  we  expect  the 


coefficients  in  the  expansion  (3.50)  to  decay  rapidly.  In  this  section  we  derive  a  bound  for  the 


rate  of  decay  of  these  coefficients. 

Lemma  3.34.  Given  an  arbitrary  choice  of  0  <  a  <  b  <  oo,  we  consider  the  SVD  of  the 


operator  (defined  in  2.61).  Let  rff  be  the  k  +  1-th  coefficient  in  the  expansion  of  the  n  +  1- 
th  left  singular  function  vn  in  the  basis  of  Laguerre  functions  (the  functions  if1*,,  defined  in 


(2.4V). 


We  define  7  =  b/ a  and  introduce  the  notation 


( 

log  2  a 

log  2b 

&max  —  UlcLX.  1 

log  7 

5 

log  7 

(3.74) 


Then, 


s 


max  : 


>  1 


(3.75) 


and 


\6k  I  <  « 


-1 

n 


J  2 

/ySmax  / 2  _ 

V  log  7 

rySmax /2  _|_ 

In  particular,  in  the  case  £7  =  C  1  ^  (as  defined  in 

Vf 


(3-5)), 


\Vk\  <  ^ 


log  7 


1  - 


l  +  v/7 


(3.76) 


(3.77) 
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Proof.  By  (3.55)  and  (2.66), 


roc  /  rb 

Vk  =  “n'  /  (  /  e_artun(t)dt  )  d>fc(w)dw. 

Jo  \Ja 


(3.78) 


Changing  the  order  of  integration  and  using  (2.62), 


Vk=an1  I  Un(t)  ((CafiT  ($*))  (t)dt. 


(3.79) 


A  simple  calculation  using  (3.11),  (3.12)  and  (3.20)  shows  that 


n  -l 


2  log  2b 
J  log  7 


V'k=an-  I  Un(s)((C1o  (Ca  b)  )(^k))(s)ds. 

1 2  log  2a 
log  7 


(3.80) 


By  the  Cauchy-Schwarz  inequality, 


kfcl  < 


-1 
n  \ 


2  log  2b 
log  7 


/  o  log  2a 
log  7 


(Cn(s))2  d-s 


1 


^2  log  2b 

log  7 


/  2  log  2a 
log  7 


((C7o(£a!6)*)(<Dfc))2(s)dS,  (3.81) 


and  by  (3.21 ), 


\Vk\  < 


Vl°g7\ 


^  log  2  b 


1087  ((c7o(£aibr)(chA:))2(S)ds. 


^  log  2a 
log  7 


(3.82) 


We  observe  that  smax  is  the  supremum  of  |s|,  where  s  G  (21j^^)2^y).  In  other  words, 


Smax  is  the  largest  magnitude  of  the  variable  s  in  the  integration  (3.82).  It  is  easy  to  observe 


that  smax  is  no  smaller  than  1.  By  (3.16), 


|((c7o((£ai6)*)($fc)))WI<« 


7 


c/2  _  1 


c/2  +  1 


(3.83) 


For  a  given  ratio  7  =  b/a ,  it  is  easy  to  observe  that  the  length  of  the  interval  ( 2 1^2° ,  2  ) 
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in  the  integral  (3.82)  is  2.  So,  by  (3.82)  and  (3.83) 


\Vk\  <  a 


_i  y/2 

\/log7 


e/2  _  ^ 
z/2  _(_  1 


(3.84) 


In  the  case  of  the  standard  form  of  the  truncated  Laplace  transform  £7,  where  a  =  &  = 

this  interval  becomes  (—1, 1),  and  smax  =  1.  So,  for  the  standard  form  of  the  truncated 
Laplace  transform, 


i y/2 


—  ftn  /i - 

Vl°g7 


7 


1/2  _  ! 


7 


V2  +  l 


=  a, 


i \/2 


\/Iog7 


1  - 


1  +  ^7 


(3.85) 


D 


Observation  3.35.  Let  vn  be  the  n  +  1-th  left  singular  function  of  Ca ^  (defined  in  (2.61[)).  Let 
vn  be  the  n  +  1-th  left  singular  function  of  (the  operator  in  the  standard  form,  as  defined  in 


(3.5)),  where  b/a  =  7.  Let  the  vectors  ryn  and  rjn  represent  the  expansions,  defined  in  (3.50), 
of  Vn  and  v„ . 


In  the  case  of  £7,  we  have  smax  =  1,  and  it  is  easy  to  observe  that  the  bound  (3.77)  for  |r/(( 


decays  faster  than  the  bound  (3.76)  for  the  general  |  fj 
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3.6  A  remark  about  the  limit  7  — >  1 


Throughout  this  dissertation  we  have  assumed  that  the  parameter  7  =  b/a  in  (3.5)  is  strictly 
larger  than  0.  In  this  section,  we  describe  some  properties  of  the  differential  operator  Du 


(defined  in  (2.82))  and  its  eigenfunctions  vn  at  the  limit  7  — >  1.  Other  aspects  of  this  limit  are 
discussed  in  [6]. 


By  substituting  b  =  a  into  (2.82), 


Du 


(w)  =  - 


d2 


du;2  du;2 


2  d  /M  +  +  (-o4^  +  2«2)  /H- 


(3.86) 


In  particular,  substituting  a  =  b  =  1/2  into  (2.82)  yields 


Du 


,  ,  d2  2  d2  Id  2  d  /  1  2  1 ' 1  .  / 

M  “  ~d^W  d^/M  +  2  dCW  d^/M  +  (  ' "l6W  +  2  )  /M- 


(3.87) 


Theorem  3.22  provides  a  relation  between  the  operator  Du  and  the  matrix  M  (see  (3.52)).  By 


substituting  a  =  b  =  1/2  into  (3.52),  we  obtain  a  diagonal  matrix: 


Mk,k  —  ~~  k(k  +  1). 


(3.88) 


Clearly,  the  eigenvalues  of  this  matrix  are 


x*  =  -k(k  +  !), 


(3.89) 


and  the  eigenvectors  are  simply  (0,  ..,0, 1,  0, ...)  .  By  theorem  3.22,  this  means  that  for  7  =  1, 


the  n  +  1-th  eigenfunction  vn  of  the  differential  operator  Du ,  is  the  Laguerre  function  <hn  (the 


function  defined  in  (2.47)),  and  the  n  +  1-th  eigenvalue  of  Du  is  Xn  =  ~n(n  +  1).  In  other 
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words,  the  Laguerre  function  4>n  is  the  solution  of  the  differential  equation 


d2  0  d2  T  Id  2  d 


1  .2  1 


"d^d^*"  +  2 d^d^  +  (  -^^  +  2  +  n(n  +  1}  )  (3'90) 


3.7  A  relation  between  the  n+  1-th  and  m+  1-th  singular  func¬ 
tions,  and  the  ratio  an/am 


Lemma  3.36.  Let  un  and  um  be  right  singular  functions,  and  let  an  and  am  be  the  correspond¬ 


ing  singular  values  of  Ca,b  (defined  in  (2.61).  Let  ifn  and  ipm  be  the  corresponding  functions 


defined  in  (2.75).  Then: 


C^n  =  Ip  ^n(x)^m(x)dx 

a n  fd  lfn(.x)lp'm(x) dx 


(3.91) 


and 


rfn  =  fa  u'n{t)um{t)dt 
01 n  fid  Wm(t)un{t)dt 


if  the  integrals  are  not  0. 


(3.92) 


Proof.  We  recall  from  (2.67)  that 


Un(t)  =  —{C(vn))(t)  =  —  f 

JO 


=  —  /  e  utvn(oj)du 


(3.93) 


Therefore,  the  derivative  of  un(t )  is 


1  /‘00 

u'n(t)  =  —  /  {—Lo)e~ultvn(u:)duj. 
JO 


(3.94) 


We  multiply  both  sides  of  the  expression  by  um(t),  integrate  both  sides,  and  change  the  order 
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of  integration: 


f  u'n(t)um(f)dt  =  —  [  (  [  ( -u)e  “^(u^dch  um(t)dt.  (3.95) 

a  ^ n  J a  0  / 


By  rearranging  the  result,  we  obtain 


u'n(t)um(t)dt  =  —  I  (- ui)vn(u)vm(uj)duj . 

JO 


(3.96) 


m  and  n  are  clearly  interchangeable,  so  that 


f 


a . 


(-uj)vn(u})vm(u)du  =  —  /  um(t)un(t)dt. 


OLr. 


(3.97) 


By  substituting  (3.97)  into  (3.96),  we  obtain  (3.92).  The  identity  (2.75)  is  used  to  obtain  (3.91). 


□ 


A  similar  relation  exists  for  the  left  singular  functions  and  their  derivatives: 


Lemma  3.37.  Let  vn  and  vm  be  left  singular  functions  and  let  an  and  am  be  the  corresponding 
singular  values.  Then: 


c*L  =  Jo°°  do; 

“n  JT  ««(wK(w)dw’ 


if  the  integrals  are  not  equal  to  0. 


The  proof  is  similar  to  the  proof  of  lemma  3.36 


(3.98) 


3.8  A  relation  between  vn(0),  and  the  singular  value  an 

The  following  lemma  provides  the  relation  between  /iq  (the  first  coefficient  in  the  expansion 


(3.33)  of  un),  un(0),  and  the  corresponding  singular  value  an. 
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Lemma  3.38.  Let  vn(uo)  be  a  left  singular  functions  of  C1  (the  operator  defined  in  (3.5)).  Let 
Uk  be  the  corresponding  right  singular  function,  and  let  an  be  the  corresponding  singular  value. 


Let  ifn  be  as  defined  in  (2.75)  and  let  hn  be  the  vector  of  coefficients  defined  in  (3.33) 
Then, 


Q.n.  — 


1 7-1  h 


2v/7  MO) 


(3.99) 


Proof.  By  the  definition  of  the  SVD  (2.66), 


{C(un))  (w)  =  anvn{Lo). 


(3.100) 


In  particular,  at  u  =  0: 


OnM0)  =  (£(lin))  (0)  =  /  un(t)dt 


(3.101) 


Using  the  change  of  variables  (2.74),  and  substituting  (2.79)  into  the  last  expression,  we  obtain: 


«n^n(0)  =  (b  —  a)  /  un(a  +  (b  —  a)x) dx  =  \J (b  —  a)  /  ipn(x)dx  (3.102) 


/  o 


'0 


Expressing  ifn  using  the  expansion  defined  in  (3.33): 


0)  =  \/(b-a)  J0  yYl  hmPm(x)j  dx 


(3.103) 


By  (2.17),  Pq{x)  =  1,  and  since  all  the  other  polynomials  Pf.  are  orthogonal  to  it, 


anvn(0)  =  V(b~  a)h 


(3.104) 
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Substituting  the  values  of  a  and  b,  defined  in  (3.3)  into  (3.104),  we  obtain  (3.99). 


Q 


3.9  A  closed  form  approximation  of  the  eigenvalues  Xn,X*nXr 
and  singular  values  an 


The  eigenvalues  of  differential  operators  Dt,  Dx  and  D,  as  the  operators  are  defined  in  (2.72), 


(2.80)  and  (2.82),  and  as  they  appear  in  equations  (2.73)  ,  (2.81)  and  (2.83),  have  closed  form 


asymptotic  expressions;  in  the  case  of  the  standard  form  £7  (the  operator  defined  in  (3.5)),  the 
eigenvalues  of  the  differential  operators  are: 


4y  _  * 

X™  _  -Q2  ~  Xn  ~  Xn  ~ 


* +  «  +  -o 

_ 

I67 


(l  +  0(n  2))  , 


(3.105) 


where  K(m )  is  the  complete  elliptic  integral  of  the  first  kind  (as  defined  in  (2.49))  .  These 


eigenvalues  are  negative  and  roughly  proportional  to  — n(n  +  1).  The  proof  for  this  asymptotic 
expression  is  involved,  and  it  will  be  provided  at  a  later  date. 

The  singular  values  an  also  have  a  closed  form  asymptotic  expression;  in  the  case  of  the 
standard  form  £7,  the  singular  values  are: 


= 


exp 


V3- 7(7 +  8x^  +  2 )K  (^) 


v/2(7  +  1) 


(l  +  0(n-1)),  (3.106) 


where  K(m )  is  the  complete  elliptic  integral  of  the  first  kind  (as  defined  in  (2.49)).  The  proof 
for  this  asymptotic  expression  is  involved,  and  it  will  be  provided  at  a  later  date. 
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Chapter  4 


Algorithms 


4.1  Evaluation  of  the  right  singular  functions  un 

In  this  section  we  introduce  an  algorithm  for  the  numerical  evaluation  of  un(t),  the  n  +  1-th 


right  singular  function  of  £7  (the  operator  defined  in  (3.5)). 


We  recall  that  un{t)  can  be  easily  calculated  from  the  function  i/jn(x)  using  (2.75).  We 


also  recall  that  ij}n{x)  is  efficiently  represented  in  the  basis  of  P*k  (the  polynomials  defined  in 


(2.16))  and  that  the  expansion  of  ipn(x)  in  P */,  is  related  to  the  n  +  1-th  eigenvector  of  the 


5-diagonal  matrix  specified  in  theorem  3.16 


The  algorithm  for  obtaining  the  right  singular  function  un(t)  is  therefore: 


Compute  hn,  the  n  +  1-th  eigenvector  of  the  matrix  M,  defined  in  (3.35). 


Compute  the  function  ipn(x)  from  hn,  using  the  expansion  specified  in  (3.33). 


Obtain  un(t)  from  'ipn(x)  using  (2.75). 


The  calculation  of  the  eigenvalues  and  eigenvectors  is  done  using  the  Sturm  sequence  method 


and  the  inverse  power  method,  as  described  in  section  2.6 
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4.2  Evaluation  of  the  left  singular  functions  vn 


In  this  section  we  introduce  an  algorithm  for  the  numerical  evaluation  of  vn,  the  left  singular 


functions  of  £7  (the  operator  defined  in  (3.5)). 


We  recall  that  vn{oj),  is  efficiently  expressed  in  the  basis  of  Laguerre  functions  as  specified 


in  (3.50).  We  recall  that  the  coefficients  of  even-numbered  left  singular  functions  of  £7  are 


given  by  the  eigenvectors  of  the  matrix  Meven,  specified  in  (3.61).  Therefore,  the  algorithm  for 


computing  an  even-numbered  left  singular  function  vn  where  n  =  2 j  is: 


Compute  r]even,J ,  the  j  +  1-th  eigenvector  of  the  matrix  Meven  specified  in  (3.59) 


Compute  the  function  v?ven(uj)  from  •qeven’:) ,  using  the  expansion  specified  in  (3.63). 


By  (3.65),  vn(uj)  =  vfen(uj) 


Similarly,  the  algorithm  for  computing  an  odd-numbered  left  singular  function  vn  where  n  = 
2  j  +  1  is: 


Compute  rjodd^ .  the  j  +  1-th  eigenvector  of  the  matrix  Modd  specified  in  (3.60). 


Compute  the  function  v°dd(cj)  from  rfdd,d ,  using  the  expansion  specified  in  (3.64). 


By  (3.65),  vn{u)  =  v°dd(uj)  . 


The  calculation  of  the  eigenvalues  and  eigenvectors  is  done  using  the  Sturm  sequence  method 


and  the  inverse  power  method,  as  described  in  section  2.6 


Remark  4.1.  Clearly,  the  left  singular  functions  of  Ca,b  (the  operator  defined  in  (2.61))  can 


be  computed  directly,  using  the  matrix  M  described  in  theorem  3.22 


4.3  Evaluation  of  the  singular  values  ar 


In  this  section,  we  introduce  two  algorithms  for  computing  the  singular  value  an  of  £7  (the 


operator  defined  in  (3.5)). 
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4.3.1 


Calculating  the  singular  value  an  from  am, 

1X371 


via  lemma  3.36  or  lemma 


Lemma  |3.36  provides  a  way  to  calculate  an+i  from  a  known  an  using  the  right  singular  func¬ 
tions.  Suppose  that  we  have  the  first  singular  value  ao-  Then,  we  calculate  the  functions  tpo 


and  -01  (the  functions  as  defined  in  (2.75))  and  use  (3.91)  to  calculate  a\  from  ao-  To  obtain 
the  other  singular  values,  we  calculate  every  an+ \  from  the  previous  an,  using  V'n+i  and 

There  are  several  obvious  methods  for  evaluating  ao  via  numerical  integration;  for  example: 


a  o 


\{C*oC)  («o))(f) 


u0(t) 


(4.1) 


and 


ao  = 


(C(u0))  (u) 

v0(u) 


(4.2) 


We  use  the  relation  provided  in  lemma  3.38  to  evaluate  ao;  in  section  4.3.2  we  use  lemma  3.38 
to  calculate  arbitrary  an  directly. 


Remark  4.2.  A  similar  algorithm,  based  on  the  left  singular  functions  vn  rather  than  the  right 


singular  functions,  is  easy  to  construct  using  lemma  3.37 


Remark  4.3.  Clearly,  if  ((£*  o  C)  ( un ))  (t)  and  un{t)  are  available  at  sufficient  precision,  re¬ 


lations  like  (4.1)  can  be  used  to  evaluate  arbitrary  singular  values.  However,  in  general,  the 


condition  number  of  the  problem  does  not  allow  high  precision  calculation  of  small  singular 


values  using  direct  integration;  and  since  an  decays  exponentially  (see  (3.106)),  only  very  few 


singular  values  can  be  calculated  by  direct  numerical  integration.  The  method  in  section  4.3.2 
provides  an  alternative  way  of  evaluating  the  singular  value  an. 
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4.3.2  Calculating  the  singular  value  an  via  lemma  (3.38) 


Let  vn(uj)  be  the  n  +  1-th  left  singular  function  of  £7  (the  operator  defined  in  (|3.5[) ) .  Let  an 
be  the  n  +  1-th  singular  value.  Let  hn  be  the  n  +  1-th  eigenvector  of  the  matrix  M  (the  matrix 


defined  in  (3.35)).  By  lemma  3.38 


Q.n.  — 


/t-  i  K 
2^7  MO)' 


(4.3) 


The  values  of  h q  and  un(0)  are  obtained  through  the  evaluation  of  the  right  and  left  singular 


functions,  as  described  in  sections  4.1  and  14.2 


Remark  4.4.  By  (3.63),  (3.64),  (2.43)  and  (2.47),  V2j(0)  is  simply  the  sum  of  the  entries  in 


the  eigenvector  r/et,enj  of  Meven  (the  matrix  defined  in  (3.59))  and  V2j+i(0)  is  simply  the  sum 


of  the  entries  in  the  eigenvector  r]odd J  of  Modd  (the  matrix  defined  in  (3.60)); 


wj<o) = 


(4.4) 


fc=0 


^•+i(o)  =  ^2% 


odd,j 


(4.5) 


k= 0 


Remark  4.5.  It  has  been  shown  in  m  that  in  some  band  matrices,  such  as  the  matrix  M  in 


(3.35),  the  first  element  of  the  vector  hn  can  be  computed  to  relative  precision,  and  not  just  to 


absolute  precision.  The  analysis  is  somewhat  involved,  and  it  will  be  reported  at  a  later  date. 
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Chapter  5 


Implementation  and  numerical 
results 


Algorithms  for  the  evaluation  of  the  right  singular  functions  un,  left  singular  functions  vn  and 
singular  values  an  of  £7  were  implemented  in  FORTRAN  77.  In  this  section,  we  present 
examples  of  numerical  experiments.  The  gfortran  compiler,  and  double  precision  arithmetic 
were  used  in  all  the  experiments,  except  for  the  last  experiment,  where  the  Fujitsu  compiler 
and  quadruple  precision  were  used. 

In  figure  [All  we  present  examples  of  right  singular  functions  un  and  left  singular  functions  vn 


of  £7  (the  operator  defined  in  (3.5)),  with  the  parameter  7  =  1.1.  The  right  singular  functions 
are  plotted  on  the  interval  (9^77  The  left  singular  functions  are  plotted  on  a  subset  of 

the  interval  (0,  00). 


Figure [572] is  the  same  as  figure  5.1,  but  with  the  parameter  7  =  10.  Figure  5.3  is  the  same 
as  figure  5.1  with  7  =  105,  and  a  different  selection  of  n. 


The  singular  values  an  of  £7,  over  a  range  of  n  and  a  range  of  7,  are  presented  in  table  5.1 


and  figure  5.4 


The  eigenvalues  x*  of  the  differential  operator  Du  (as  defined  in  (2.82))  are  presented  in 
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table  |5.2|  and  figure  |5.5[ 

we  plot  of  un{ a);  the  right  singular  function  of  £7,  evaluated  at  the  point 
we  plot  un(0);  the  left  singular  function,  evaluated  at  the  point  u  =  0. 


In  figure 


5.6 


5.7 


a  =  In  figure 

An  analysis  of  the  properties  of  un  and  vn  at  these  endpoints  will  be  presented  at  a  later  date. 
In  table  5.3  we  present  several  singular  values  smaller  than  lO-1000. 
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(b)  Left  singular  functions  vn. 

Figure  5.1:  Singular  functions  of  £7,  where  7  =  1.1. 
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(a)  Right  singular  functions  un. 
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(a)  Right  singular  functions  un. 


(b)  Left  singular  functions  vn. 

Figure  5.2:  Singular  functions  of  £7,  where  7  =  10. 
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04 


n  = 


(a)  Right  singular  functions  un. 


(b)  Left  singular  functions  vn. 

Figure  5.3:  Singular  functions  of  £7,  where  7  =  105. 
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Figure  5.4:  Singular  values  an  of  C ~ 
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Figure  5.6:  un  (l/(2y^y)).  The  right  singular  functions,  evaluated  at  t  =  a  =  1/(2^). 
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Figure  5.7:  un(0).  The  left  singular  functions,  evaluated  at  uj  =  0  . 
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Table  5.1:  Singular  values  an  of  C, 
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Table  5.2:  Eigenvalues  %*  of  the  differential  operator  Du  (defined  in  2.82). 
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Table  5.3:  Examples  of  singular  values  an  smaller  than  10 


7 

n 

Qt-n 

1.1E  +  0 

520 

8.70727 E  -  1002 

1.0E  +  1 

1721 

3.6693415  -  1001 

1.0E  +  2 

2797 

5.2996115  —  1001 

1.0E  +  3 

3872 

5.7114 6E  -  1001 

1.0E  +  4 

4946 

9.4419115  -  1001 

1.0E  +  5 

6021 

8.8974815  -  1001 

1000 
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Chapter  6 


Conclusions  and  generalizations 


In  this  dissertation  we  have  introduced  efficient  algorithms  for  the  evaluation  of  the  singular 
functions  and  singular  values  of  the  truncated  Laplace  transform. 

Among  the  obvious  generalizations  of  this  work,  is  the  Laplace  transform  in  higher  dimen¬ 
sions.  Another  closely  related  object  is  the  two-sided  band-limited  Laplace  transform,  £c;  for 
a  given  c  6  C  and  a  function  /  6  L2(— 1, 1),  the  later  is  defined  by  the  formula 

(4(/))  M  =  J  e~ct“f(t)dt.  (6.1) 

As  we  will  report  in  more  detail  at  a  later  date,  much  of  the  analysis  of  Tc  (the  operator  defined 
has  a  natural  extension  to  Cc. 

One  of  the  results  of  this  work  will  be  the  construction  of  interpolation  formulas  in  the  span 
of  right  or  left  singular  functions,  as  well  as  associated  quadrature  formulas. 

In  a  future  paper  we  will  discuss  asymptotic  properties  of  the  truncated  Laplace  transform 
and  of  the  associated  differential  operators,  and  the  relations  between  all  these  operators. 


in  (1.3 
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